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FUSE OBSERVATIONS OF INTRINSIC ABSORPTION IN THE SEYFERT 1 GALAXY MRK 509 
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ABSTRACT 

We present far-ultraviolet spectra of the Seyfert 1 galaxy Mrk 509 obtained in 1999 November with the Far 
Ultraviolet Spectroscopic Explorer (FUSE). Our data span the observed wavelength range 915-1 185 A at a res- 
olution of 20 kms~'. The spectrum shows a blue continuum, broad O VI AA1032, 1038 emission, and a broad 
C III A977 emission line. Superposed on these emission components, we resolve associated absorption lines of 
O VI AA1032, 1038, C III A977, and Lyman lines through L^. Seven distinct kinematic components are present, 
spanning a velocity range of -440 to + 170 kms~' relative to the systemic velocity. The absorption is clustered 
in two groups, one centered at -370 kms"' and another at the systemic velocity. The blue-shifted cluster may be 
associated with the extended line emission visible in deep images of Mrk 509 obtained by Phillips et al. Although 
several components appear to be saturated, they are not black at their centers. Partial covering or scattering per- 
mits ~ 7% of the broad-line or continuum flux to be unaffected by absorption. Of the multiple components, only 
one has the same ionization state and column density as highly ionized gas that produces the O VII and O VIII 
ionization edges in X-ray spectra of Mrk 509. 

Subject headings: Galaxies: Active — Galaxies: Individual (Mrk 509) — Galaxies: Nuclei — Galaxies: 

Quasars: Absorption Lines — Galaxies: Seyfert — Ultraviolet: Galaxies — X-Rays: Galaxies 



1. INTRODUCTION 

Mrk 509 straddles the boundary in luminosity between 
Seyfert 1 nuclei and QSOs with an absolute magnitude of 
Mb = -22.0. This gives it particular importance for understand- 
ing how the properties typical of the well-studied nearby active 
galactic nuclei (AGN) might scale with luminosity, and, there- 
fore, be applied to the more distant and luminous QSOs. Like 



tion Unes (York et al. 1984; 


Crenshaw, Boggess, & Wu 1995 


Savage, Sembach, & Lu 1997 


; Crenshaw et al. 19991) and X-ray 



et al. 1998). The UV absorption is blue-shifted relative to the 
systemic velocity of 10,365 kms"' ( [Phillips et al. 1983), and 



Phillips et al. suggest that it is related to the expanding shell of 
high ionization gas visible in images and spectra out to a radius 
of 15" from the nucleus. 

At high resolution, UV absorption in Seyferts appears to 
be kinematically complex (NGC 3516: Crenshaw, Maran, & 
Mushotzky 1998; NGC 5548: [Mathur, Elvis, & Wilkes 1999| ), 
with at most one of the UV components having any possible 
relati on to the X-ray absorbing gas (M athur, Elvis, & Wilkes 
1999). High resolution observations of UV absorption lines 
covering a range of ionization states is the key to determin- 
ing which, if any, of the UV absorbers are associated with an 
X-ray warm absorber. The FUSE observations of Mrk 509 
we present in this paper extend the far-UV spectral coverage 
to wavelengths shortward of 1200 A. Our spectra include the 
O VI AA1032, 1038 resonance doublet and the high-order Ly- 



man lines, down to the redshifted Lyman limit. Since the 
La lines resolved in earlier UV observations are saturated, the 
high-order lines visible with FUSE provide a better constraint 
on the total neutral hydrogen column density. The O VI doublet 
is a crucial link for establishing a connection between the higher 
ionization absorption edges seen in the X-ray and the lower ion- 
ization absorption lines seen in earlier UV observations. Based 
on these new observations, we discuss the implications for the 
location of the absorbing gas in AGN, and how intrinsic absorp- 
tion in low-z AGN may be related to the broad-absorption line 
phenomenon in more luminous QSOs. 

2. OBSERVATIONS 

FUSE comprises four separate primary mirrors gathering 
light for four prime-focus, Rowland-circle spectrographs and 
two, two-dimensional, photon-counting detectors. For a full 
description of FUSE, its mission, and its in-flight performance, 
see Moos et al. (2000), and Sahnow et al. (2000). Two of 
the optical systems employ LiF coatings on the optics, giving 
coverage from ^ 990-1 187 A, and the other two use SiC coat- 
ings, which provide reflectivity down to wavelengths as short 
as 905 A. The mirror systems focus light on a slit assembly in 
the focal plane. The holographically ruled gratings disperse the 
light entering through the slits and form an astigmatic image 
on the two-dimensional microchannel-plate detectors. The de- 
tectors have KBr photocathodes and delay-line anode readouts 
that provide the location and arrival time of each photon event. 

Mrk 509 was observed on 1999 Nov 9 and 1999 Nov 11 
through the 30" x 30" low-resolution apertures. We obtained 
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Fig. 1 . — FUSE spectrum of Mrk 509 (binned by 20 pixels) is shown as a thin black line. Hopkins Ultraviolet Telescope (HUT) data obtained during the Astro-2 
mission in 1995 March (scaled by a factor of 0.69) are shown as the green curve. 



good spectra from the LiFl and LiF2 channels covering the 
987-1187 A band, and lower signal-to-noise-ratio (S/N) spec- 
tra from the SiC2 channel covering 905-1091 A. As a result 
of channel misalignment during the observations, almost no 
data were obtained in the SiCl channel. We recorded the data 
in photon-address mode, which provides a time-tagged list of 
event positions in the down-linked data str eam. This enables us 
to filter out short periods of event "bursts" (Sahnow et al. 200C ) 
and to correct for image motion on the detectors to achieve 
the best spectral resolution. As described by Sahnow et al. 
(2000), we extracted one-dimensional spectra from the two- 
dimensional data on each of the active detector segments. These 
extracted spectra are dark subtracted and flux and wavelength 
calibrated by the standard FUSE calibration pipeline. We es- 
timate that the flux scale is accurate to ^10%, and that wave- 
lengths are accurate to ^15 kms"'. Poisson errors and data 
quality flags are propagated through the data reduction process 
along with the science data. 

To produce the full spectrum shown in Figure we have 
spliced together sections of data from segments SiC2A, LiFl A, 
SiC2B, LiFlB, and LiF2A to eliminate the gaps in wavelength 
coverage and to present the best S/N. These data have also been 
binned by 20 pixels (0. 12 A) to show the overall appearance of 
the spectrum. In our analysis discussed in the following sec- 
tions, we use spectra binned by only 5 pixels. This preserves 
the full spectral resolution of ^ 20 kms~' for this observation. 

In Figure |l] strong, broad O VI emission is readily apparent 
as is broad C HI A977 emission. Another less prominent fea- 
ture in the spectrum is the hump of emission redward of O VI 
AA1032, 1038. We have marked the position of He II A1085, 
but one can see that its wavelength is a bit too long. A simi- 
lar feature was noted in the spectra of low-redshift quasars by 
Laor et al. (1995) who suggested this may be blended Fe II 
emission. The numerous Galactic absorption features, particu- 
larly H2, make the intrinsic spectrum difficult to see shortward 
of 1000 A, but, at the resolution of FUSE, one can trace out 
the intrinsic spectrum as the peaks between the foreground ab- 
sorptions. Note how the C III A977 line is not prominent in the 
overlayed Hopkins Ultraviolet Telescope (HUT) data, primar- 
ily because the foreground Galactic absorption renders it less 
visible at the ^ 3 A resolution of HUT. 



3. WARM ABSORBING GAS IN MRK 509 

A detailed examination of the FUSE spectrum at full spectral 
resolution shows that absorption near the redshift of Mrk 509 
is visible in the Lyman lines, in the O VI AA1032,1038 res- 
onance doublet, and in C III A977. Figure 2 shows portions 
of the Mrk 509 spectrum around each absorption line complex 
with velocities relative to the AGN systemic velocity. 

Close inspection of the O VI doublet and the Lyman lines 
shows that while much of the absorption appears to be satu- 
rated, the absorption troughs are not bla ck. (Scattered light i n 
FUSE at these levels is negligible. See Sahnow et al. 2000 ) 
Partial covering by the absorbers or scattering around the ab- 
sorbing region can explain this appearance. Extended broad- 
line Balmer emission and varia ble, polarized broad Ba lmer-line 
emission are seen i n Mrk 509 (Mediavilla et al. 1998 ^Thomp- 
son & Martin 1988; Young et al. 199'y ). These are both indica- 



tions that scattering may play some role in producing the light 
at the bottoms of the absorption troughs. 

In the L/3 and O VI absorption lines we have identified at least 
7 distinct kinematic components spanning a range of velocities 
from -478 to +166 kms"' relative to the systemic velocity of 
Mrk 509. To ascertain the physical properties of these intrinsic 
absorbers in Mrk 509, we fit a model to the spectral regions sur- 
rounding the O VI lines, t he Lyman li nes, and C III A977 using 
the IRAF task specf it ( Kriss 1994 ). Figure 3 shows the best 
fit overlayed on the L/3/0 VI region. Our model of the emis- 
sion components includes an underlying power-law continuum, 
a pair of broad (FWHM'^l 1,000 kms"') O VI emission lines 
with their relative intensities fixed at the optically thin ratio of 
2:1, a pair of narrow (FWHM^1200 kms"') O VI emission 
lines, again with their relative intensities fixed at a 2:1 ratio, 
and a broad L/3 emission line with its width and velocity linked 
to those of the broad O VI lines. The absorption lines are treated 
as Gaussians in optical depth, and they are allowed to partially 
cover the emission components with fraction fc. The wave- 
lengths of the O VI doublets are linked at the ratio of their labo- 
ratory values, their velocity widths are required to be identical, 
and their relative optical depths are fixed at a 2:1 ratio. Thus, 
a pair of lines uniquely determines the column density and the 
covering fraction for a given kinematic component. The L/3 
lines have their wavelengths linked to those of the O VI lines, 
but, to allow for residual uncertainties in the FUSE wavelength 
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scale, we permit a linear adjustment to the whole group of L/3 
lines. The widths and the covering fractions for the L/3 lines are 
fixed at the values determined for the O VI lines. 
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Fig. 2.— Normalized line profiles for the Lyman lines, C III A977, O VI A1032 
(OVIb), and O VI A1038 (OVIr) ai'e shown. These data are binned by 5 pixels, 
and foreground Galactic absorption features have been divided out. The central 
wavelengths of the 7 distinct intrinsic absorption components are marked. The 
La profile at the upper left is archival HST data from a GHRS observation by 
B. Savage. 

Higher-order Lyman lines are visible in the FUSE spectrum 
out to L^. These lines are also fit simultaneously with L/3. Their 
optical depths are fixed at the ratio of their oscillator strengths 
to that of L/3, their wavelengths are linked to those of O VI at 
the ratio of the vacuum wavelengths (again allowing for slight 
linear adjustments), and their widths are fixed at the same val- 
ues as those of the O VI components. The Lyman lines permit 
an independent check of the covering fractions determined from 
the O VI doublets. We find that if they are allowed to vary in- 
dependently, they give results consistent with those determined 
from O VI alone, so we leave them fixed at the O VI values in 
our final results. C III A977 is treated similarly. Although there 
is no independent check on its covering fraction, we also fixed 
/c at the O VI values. 

The results of our fits for C III A977, L/3, and O VI A 1032 
are given in Table 1. The column densities for each ion are 
determined by integrating the optical depth across each param- 
eterized line profile. (Since the parameters given completely 
determine the values for O VI A1038 and the remaining Lyman 
lines, these are not shown.) 

To determine physical conditions in the absorption compo- 
nents, we used photoionization models similar to those used by 
Krolik & Ki-iss (1995) and Ki'iss et al. (1996). From a grid 
of models we determined the total column density and the ion- 
ization parameter based on the observed relative columns of 
H I and O VI. With no other constraints, this method can lead 
to double-valued results for the ionization parameter and col- 
umn density since the ratio of O VI to H I will rise to a peak 
and then decline. However, we note that the presence of C III 
in components 1-4 restricts the solutions to the lower ioniza- 
tion parameter in these cases. Component 5 lies near the peak 
of the N(0 Vl)/N(H l) curve, and, as we show below, the in- 
ferred ionization and column density are corroborated by X-ray 



observations of O VII and O VIII. As component 5 fully ac- 
counts for the observed X-ray absorption, a high-column, high- 
ionization-parameter solution for components 6 and 7 is also 
ruled out by the X-ray observations. Physical parameters for the 
7 kinematic components are shown in Table 2. Note that most 
components have relatively low total column densities and ion- 
ization parameters. The components associated with the most 
blue-shifted complex could well be associated with the outflow- 
ing narrow-emission line gas as originally suggested by Phillips 
et al. (1983). This would place it many kiloparsecs from the 
central ionizing source. Similarly, the lower ionization com- 
ponents (6 & 7) near the systemic velocity may be associated 
with the low -ionization gas in th e rotating disk near the center 
of Mrk 509 ( ^hillips et al. 1983| ). 
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Fig. 3.— FUSE spectrum of Mrk 509 (binned by 5 pixels) in the L/3/0 VI 
region is shown as the thin black line. The best fit described in the text is over- 
layed in color. The thin green line shows the fitted continuum and emission 
components. The thin red line shows the fitted intrinsic absorption components, 
and the thin blue line shows the foreground Galactic absorption lines. The cen- 
tral wavelengths of the 7 distinct intrinsic absorption components are marked. 

Among all the absorption components, the most exceptional 
is #5. Its high ionization parameter and high total column den- 
sity make it likely to be directly associated with the intrinsic 
X-ray warm absorber In his analysis of the ASCA X-ray spec- 
trum of Mrk 509, Reynolds (1997) found O VII and O VIII op- 
tical depths of 0.11:l:o o| and 0.04+q q3, respectively. For thresh- 
old photoionization cross secti ons of gp? = 0.239 x 10"' ^ cm^ 
and fTo8 = 0.109 x 10"'** cm- ( ^eilman & Manson 1979[ ), this 
implies column densities of Noi = (4.6;!;j 7) ^ 1*^" cm and 
■^08 = (3.7;!;2g) X 10'^ cm"^. Our photoionization modeling for 
UV absorption component 5 predicts Noi = 2.3x10'^ cm"^ and 
Noa = 0.1 X 10'^ cm~^. Considering the large uncertainties in 
the X-ray columns and the temporal difference between the X- 
ray and UV observations, the agreement is remarkably close. 
The X-ray absorption predicted by our photoionization model- 
ing for the other UV absorption components is negligible. 

We conclude that observations of the O VI doublet in AGN 
is a perfect complement to observations of the O VII and O VIII 
edges often seen in the X-ray spectra of Seyfert 1 galaxies. 
Comparison of the UV and X-ray measurements permits us to 
rigorously test whether the gas responsible for X-ray absorption 
in Seyferts also gives rise to the UV absorption, as suggested 
by Mathur et al. (1994, 1995). The multiple kinematic compo- 
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nents, the wide range of ionization parameters, and the typical 
low total column s inferred for most of t he UV absorbers (see 
Kriss et al. 1996 ; Crenshaw et al. 1999 ) make it unlikely that 



the UV and X-ray absorption arise in the same gas. In fact, in 
Mrk 509, with the high spectral resolution and far-UV sensitiv- 
ity of FUSE, we can easily pick out the high ionization absorp- 
tion component that is likely to be directly associated with the 
X-ray absorbing gas. The remaining components (which dom- 
inate the total UV absorption) are lower column density and 
lower ionization. This suggests that the absorbing medium is 
complex, with separate UV and X-ray dominant zones. One 
potential geometry is high density, low column UV-absorbing 
clouds embedded in a low density, high ionization medium that 
dominates the X-ray absorption. This is possibly a wind driven 
off the obscuring torus or the accretion disk. 

TABLE 1 
Absorption Lines in Mrk 509 



Feature 


# 


(A) 


Nion 

(cm-2) 


Av" 
(kms-') 


FWHM 

(kms-') 


fc 


cm 


1 


0.11 


1.6 X 10" 


^38 


49 


0.77 


A977.02 


2 


0.21 


2.4 X 1013 


-349 


63 


0.91 




3 


0.06 


5.4 X lO'^ 


-280 


31 


1.00 




4 


0.17 


2.1 X 10" 


-75 


57 


0.92 




5 


<0.01 


< 2.0 X 10'2 


-5 


40 


0.92 




6 


<0.01 


< 2.0 X 10'2 


+71 


48 


0.94 




7 


< 0.01 


< 2.0 X 10'2 


+166 


59 


0.94 


HI 


1 


0.25 


5.7 X lO'* 


-438 


49 


0.77 


A1025.72 


2 


0.24 


5.3 X lO'"* 


-349 


63 


0.91 




3 


0.26 


9.3 X 10'* 


-280 


31 


1.00 




4 


0.43 


6.0 X 10''^ 


-75 


57 


0.92 




5 


0.36 


2.7 X 10'* 


-5 


40 


0.92 




6 


0.25 


1.2 X 10'5 


+71 


48 


0.94 




7 


0.28 


1.2 X 10'5 


+166 


59 


0.94 


O VI 


1 


0.24 


1.9 X lO'* 


-438 


49 


0.11 


A1031.93 


2 


0.46 


1.3 X lO''' 


-349 


63 


0.91 




3 


0.14 


1.2 X 10'* 


-280 


31 


1.00 




4 


0.44 


1.2 X 10'5 


-75 


57 


0.92 




5 


0.31 


3.2 X 10'5 


-5 


40 


0.92 




6 


0.37 


1.2 X 10'^ 


+71 


48 


0.94 




7 


0.41 


9.4 X 10'* 


+166 


59 


0.94 



"Velocity relative to a systemic redshift of cz - 10365 km s ' (Phillips et al. 
1983). 

One puzzling inconsistency with the outflow hypothesis, 
however, is that component 5 is at rest with respect to the 
systemic velocity. If the X-ray absorbing gas is in the low- 
ionization disk observed by Phillips et al. (1983), then the low 
line-of-sight velocity could be explained by having most of its 
motion transverse to the line of sight. The ionization parame- 



ter U ^ 0.4 is consistent with this location if the gas density is 
low enoug h— for n = 10^ cm ^ and an ionizing luminosity of 
L,„„ = 3.4 X 10'*'' erg s"' for Mrk 509, the absorbing gas would 
be located at a distance of 330 pc from the nuclear source, or 
~ 0.5" for H() = 65 km s"' Mpc"'. Setting better constraints on 
the location of the X-ray absorbing gas, however, will require 
better knowledge of the gas density, which can be determined 
from variability studies ( e.g., [Kriss et al. 1997|; H amann, Bar- 
low, & Junkkarinen 1997; |Espey et al. 1998| ). 

TABLE 2 

Physical Properties of the Absorbers in Mrk 509 



# 


Novi/Nhi 




logU 






(cm--) 




1 


0.37 


2.0 X 10"* 


-1.64 


2 


1.51 


5.9 X 10"* 


-1.19 


3 


0.48 


2.2 X 10"* 


-1.79 


4 


0.19 


1.6 X 10" 


-1.73 


5 


13.9 


5.0 X lO^o 


-0.43 


6 


2.14 


7.6 X 10"* 


-1.41 


7 


2.76 


6.8 X 10"* 


-1.46 



4. SUMMARY 

The far-UV spectrum of the Seyfert 1 galaxy Mrk 509 shows 
bright broad O VI AA1032,1038 emission as well as broad C III 
A977 emission. Kinematically complex intrinsic absorption in 
Mrk 509 shows at least 7 distinct components in the Lyman 
lines, O VI AA1032,1038, and Cm A977. Although many of 
the O VI and Lyman line components appear to be saturated, 
they are not black, implying that partial covering or scatter- 
ing affects the absorption. Only one of the intrinsic absorption 
components in Mrk 509 is likely to be associated with the warm 
X-ray absorbing gas. Component 5 (near the systemic velocity) 
has an ionization state and column density that is in reasonable 
agreement with the O VII and O VIII absorption edges seen in 
the ASCA X-ray spectrum. The high resolution of FUSE and 
its sensitivity in the O VI band make it an ideal tool for identify- 
ing high-ionization UV absorbers that may correspond to X-ray 
warm absorbers. 

This work is based on data obtained for the Guaranteed Time 
Team by the NASA-CNES-CSA FUSE mission operated by the 
Johns Hopkins University. Financial support to U. S. partici- 
pants has been provided by NASA contract NAS5-32985. G. 
Kriss acknowledges additional support from NASA Long Term 
Space Astrophysics grant NAGW-4443. 
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